Observation of exotic meson production in the reaction tt p — > t/'-k p at 18 GeV/c 
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An amplitude analysis of an exclusive sample of 5765 events from the reaction n~p — » r]'n~p at 
18 GeV/c is described. The n'ir~ production is dominated by natural parity exchange and by three 
partial waves: those with J PC = + , 2 ++ , and 4 ++ . A mass-dependent analysis of the partial-wave 
amplitudes indicates the production of the 02 (1320) meson as well as the 04 (2040) meson, observed 
for the first time decaying to n'-K~ . The dominant, exotic (non-gq) 1 v partial wave is shown to be 
resonant with a mass of 1.597 ± O.OlOl^ow GeV/c 2 and a width of 0.340 ± 0.040 ± 0.050 GeV/c 2 . 
This exotic state, the 7Ti(1600), is produced with a t dependence which is different from that of the 
02(1320) meson, indicating differences between the production mechanisms for the two states. 



Exotic mesons - those whose valence structure is not 
composed of a quark-antiquark (qq) pair - have been dis- 
cussed Jl 10 for many years but have only recently been 
observed experimentally. The underlying structure of the 
observed exotic states at 1.4 GeV/c 2 decaying into r]ir~ 
0-|l3| and at 1.6 GeV/c 2 decaying into p°ir- JlJ] is not 
yet understood. Possible explanations for these 1=1 
states could be that they are hybrid mesons, consisting 
of a qq pair and a constituent gluon, or four-quark (qqqq) 
states. However, within the framework of the flux-tube 
model the masses of these states are somewhat low to be 
hybrid mesons Jq] ; and four-quark states are expected to 
be very broad (lj. 

Since the models for exotic mesons typically predict 
masses, widths, and branching ratios, and since it is im- 
portant to classify the exotic states to provide necessary 
input to the QCD models, there is a strong motivation to 
search for additional states as well as to search for addi- 
tional decay modes for the observed states. In this paper, 
we describe the search for exotic states decaying into the 
r/'n^ final state using the reaction ir~p — > r/'n^p, where 
rj' — > ?77r + 7r~ and r\ — » 77. 

The data sample was collected during the 1995 run 
of experiment E852 at the Multi-Particle Spectrometer 
facility at Brookhaven National Laboratory (BNL). A 7r~ 
beam with laboratory momentum 18 GeV/c and a liquid 
hydrogen target were used. A detailed description of the 
experimental apparatus can be found elsewhere |H| . 

The trigger required three forward-going charged 
tracks, a charged recoil track and a signal in a lead-glass 
electromagnetic calorimeter (LGD). A total of 165 mil- 
lion triggers of this type were recorded. After reconstruc- 
tion, 1.37 million events satisfied the trigger topology and 
had two clusters in the LGD. The r\ signal is seen in 
the 77 effective mass distribution in Figjj] (a). Apply- 
ing kinematic fitting some 70000 events consistent 
with the jM77r + 7r~7r~ (77 — > 77) final state were found. 
These events satisfied energy-momentum conservation at 
the production and 77 decay vertices with a confidence 
level CL > 0.05 as well as the requirement that the dif- 
ference between the azimuthal angles of the the fitted 
proton direction and the measured proton track be less 
than 10°. As seen from the r)n + n~ effective mass dis- 
tribution (uncorrected for acceptance) in Fig.|] (a) the rj' 
signal lies over an approximately 10% non-77' background. 
The second peak in the 7?7r + 7r~ mass spectrum is due to 
production of the /i(1285) and ry(1295) resonances. 

The next level of selection identified 6040 events consis- 
tent with the pr]'ir~ (rj' — > i]ir + ir~ , 77 — > 77) final state. 
These events satisfy energy-momentum conservation at 
the production, 7/ and 77 decay vertices with CL > 0.05 
as well as topological and fiducial volume cuts. The re- 
sulting uncorrected r]'Tr~ mass spectrum (Figj^ (b)) has 
a broad peak near 1.6 GeV/c 2 and structure around 1.3 
GeV/c 2 . 



momentum-transfer \t\ is shown in Fig.g (a). The am- 
plitude analysis discussed below was made for the data 
in the range 0.09 < \t\ < 2.5 GeV 2 /c 2 . Because of the 
very low acceptance in the region \t\ < 0.09 GeV 2 /c 2 , 
the 275 events in that region were not used. In the in- 
terval 0.25 < \t\ < 1.0 GeV 2 /c 2 the \t\ distribution has 
an exponential behavior and can be fitted with the func- 
tion f(t) = ae b l*l with 6= -2.93 ±0.11 (GeV/c)" 2 . The 
magnitude of b is significantly less than that observed for 
the r/ir~ final state [jll 12|, where b s» —5 (GcV/c)~ 2 (see 
the discussion below). 

A mass-independent partial-wave analysis (PWA) 
Jl2],[L6 17 j of the data was used to study the spin-parity 
structure of the tj'tt" system. The partial waves are pa- 
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the angular momentum, P the parity and C the C parity 
of the rj'ir~ system; m is the absolute value of the angular 
momentum projection; and e is the reflectivity (coincid- 
ing with the naturality of the exchanged particle ijLsf). 
We will use simplified notation in which each partial wave 
is denoted by a letter, indicating the rfir^ system's angu- 
lar momentum in standard spectroscopic notation, and 
a subscript, which can take the values 0, +, or — , for 
m e = Q~, 1 + , or 1~ respectively. We assume that the 
contribution from partial waves with m > 1 is small and 
can be neglected. 

PWA fits shown 
0.05 and 0.10 
2.5 GeV/c 2 and 
-, D , D + and G + 
For each partial wave the complex production amplitudes 
are determined from an extended maximum likelihood fit 
p7[ . The spin-flip and spin-nonflip contributions to the 
baryon vertex lead to a production spin-density matrix 
with maximal rank two. A rank two mass-independent 
PWA in a system of two pseudoscalars cannot be per- 
formed because of the presence of a continuous mathe- 
matical ambiguity. Rank-two fits were done when addi- 
tional assumptions for the amplitudes were introduced 
(assumptions regarding the t-dependence and mass de- 
pendence of the amplitudes) to resolve the continuous 
ambiguity problem, and they gave results consistent with 
those from the rank one fits. The PWA fits presented in 
this paper are with spin-density matrix of rank one. 

The experimental acceptance was determined by com- 
parison of the data with a Monte Carlo event sample. 
The Monte Carlo events were generated with isotropic 
angular distributions in the Gottfried-Jackson frame. 
The detector simulation was based on the E852 detector 
simulation package SAGEN |Tl|Jl^ ]. The experimental 
acceptance was incorporated into the PWA by means of 
Monte Carlo normalization integrals |i~2| . The quality of 
the fits was determined by a x 2 comparison of the exper- 
imental multipole moments with those predicted by the 
results of the PWA fit ill. 
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GeV/c 2 fits and Fig. | for the 0.10 GeV/c 2 fits. The 
former are intended to show detail in the high statistics 
low-mass region and the latter are used to study the high- 
mass region. The unnatural-parity-exchange waves (not 
shown) are small, poorly determined, and do not affect 
our results. 

The acceptance-corrected numbers of events predicted 
by the PWA fits for the stronger partial waves and the 
phase differences between them are shown as points with 
error bars. There are discrete mathematical ambiguities 
in the description of a system of two pseudoscalar mesons 
fj"7f . The ambiguous solutions were found by performing 
1000 PWA fits with random starting values in each mass 
bin. The range of the ambiguous solutions in a mass 
bin is presented by black rectangles, and the maximum 
extent of their statistical errors is shown as the error 
bar. In most mass bins, the range of values found for the 
ambiguous solutions was small enough that they cannot 
be distinguished in the figures. 

Between 1.5 and 1.8 GeV/c 2 the exotic P+ is the dom- 
inant wave. Its intensity distribution consists of a broad 
structure peaked near 1.6 GeV/c 2 . The D + wave inten- 
sity has a narrow peak at 1.3 GeV/c 2 associated with 
the 02(1320) and a broad structure at higher masses. 
The G+ wave intensity is negligible below 1.7 GeV/c 2 
and is clearly nonzero in the higher mass region. The 
(P+ - D + ), (D + - G+) and (P+ - G+) wave phase 
differences show rapid changes possibly indicative of the 
presence of interfering resonant states. The observed P+ 
and D + intensities and their relative phase difference dis- 
tribution are consistent with those reported by the VES 
collaboration. Leakage studies |l2| were carried out 
and no leakage of significance was found among the dom- 
inant waves. 

To study the nature of the observed partial waves, 
three different kinds of mass-dependent analyses (MDA) 
flijf have been carried out. In the first fit type (Fit 1), the 
P + and D + intensities and their phase difference were fit- 
ted using the PWA results in 0.05 GeV/c 2 bins. For the 
other two types (Fit 2 and Fit 3), the P+ and G+, and 
the D + and G+ waves respectively were fitted, along with 
their phase differences, using the 0.10 GeV/c 2 PWA re- 
sults. These fits all used linear combinations of relativis- 
tic Breit-Wigner functions (poles) with mass-dependent 
widths and Blatt-Weisskopf barrier factors. 

Because of the presence of distinct ambiguous PWA 
solutions in some mass bins, all possible combinations 
of these solutions were used as inputs to the MDA fits. 
Typical fits are shown as the smooth curves in Figs. || 
and [| The fits use a single Breit-Wigner function to 
describe the P+ wave, and two Breit-Wigner functions 
to describe both the D + and G+ partial waves. (The 
fits can be improved by a second P+ resonance in the 
1.4 GeV/c 2 region corresponding to the 7r(1400) state 
observed previously |Tl]-|T3|]. However, since the fits are 
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final state is not required in the analysis.) 

Many acceptable mass-dependent fits {x 2 /dof < 1.5) 
were obtained for each of the three fit types. The reso- 
nance parameters from these fits were all consistent with 
each other. Results from all of these fits were thus re- 
tained to determine resonance parameters. The reso- 
nance parameters are given in Table [j| for the P+ ex- 
otic resonance and for the lower-mass resonances in the 
D + and G+ waves. In this table, the central value for the 
mass and width of each resonance as well as the statistical 
error in these quantities were determined as the average 
of those quantities over all acceptable fits. The first error 
in these values is statistical, determined using the covari- 
ance matrix of the mass-independent PWA; the second 
is systematic. The systematic errors are based on the 
range of values allowed by taking into account different 
assumptions for the partial widths of the states, different 
parameterizations of the D + wave, and different ambigu- 
ous solutions. The experimental resolution has not been 
unfolded. 

The mass and the width of the P+ state are consistent 
with those of the 7Ti(1600) exotic state observed in the 
7r + 7r~7r~ system Ji~4| . Our data are thus consistent with 
the observation of a second decay mode of the 7Ti(1600). 

The first pole in the D + partial wave has mass and 
width consistent with those of the ci2(1320). The D + 
wave in the mass region above 1.4 GeV/c 2 is consistent 
with a broad Breit-Wigner function centered around 1.8 
- 1.9 GeV/c 2 and with width between 0.55 and 0.75 
GeV/c 2 . Alternatively, that mass region can be de- 
scribed with two narrower Breit-Wigner functions. The 
alternative parameterizations of the D + wave do not af- 
fect the conclusions of this paper. 

The G+ partial wave has been parameterized with two 
Breit-Wigner functions in fit types 2 and 3. The parame- 
ters of the lower-mass state given in Table [Hare consistent 
with the mass and width of the 04 (2040) |2lf|, which has 
not been observed previously in the r/'n" system. The 
second Breit-Wigner pole in the G+ wave is located in the 
high-mass region (« 2.4 GeV/c 2 ) where limited statistics 
results in sizable statistical and systematic uncertainty. 
Its physical interpretation is unclear. 

Mass-independent analyses have also been performed 
for two separate four-momentum-transfer intervals in- 
cluding equal numbers of events (see Fig. |J). The fits 
show that, as \t\ increases, the production rate for the 
a2(1320) decreases faster than the production rate for 
the exotic state. Note from Fig. || that the numbers of 
events in the marked peak bins for P + production are 
nearly equal for the two \t\ intervals (the ratio of the 
high-|i| to the low-|i| peaks is 1.00 ± 0.12) whereas for 
D + production, the ratio is 0.71 ± 0.15. Since the \t\ 
distribution is correlated with the production mechanism 
for peripheral processes, we conclude that exotic meson 
production proceeds via a different production mecha- 
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or that it proceeds with a different mixture of the same 
production mechanisms. 

In conciusion, we have studied the rji(~ system pro- 
duced in the reaction ir~p — > pij'ir^ at 18 GeV/c . We 
find that an exotic meson, the 7Ti(1600) is produced, de- 
caying to t/tt - . The different t dependence for their pro- 
duction shows that the well-known 02 (1320) and the ex- 
otic tti (1600) are produced via different production mech- 
anisms. Finally, a high-mass state consistent with the 
a 4 (2040) has been observed decaying to t]'tt~ . 
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TABLE I. Fitted Resonance Parameters 



Partial Wave 


Mass 


Width 


P+ 

D+ 
G+ 


1.597 ±0.010^^;^o 
1.318 ± 0.008l°:g O g 
2.000 ± 0.040lo.o2o 


0.340 ± 0.040 ± 0.050 
0.140 ±0.035 ±0.020 
0.350 ±0.100to.o5o 
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FIG. 1. (a) The r]ir + ir~ effective mass distribution for 
events consistent with the reaction ir~p — > prj-K + n + (two 
entries per event). The inset shows the 77 effective mass dis- 
tribution in 0.01 GeV/c 2 bins, (b) The Tj-K effective mass 
distribution. The distributions are uncorrected for accep- 
tance. The smooth curve in (b) shows the true mass accep- 
tance based upon the angular distributions determined in the 
partial wave analysis. 
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FIG. 2. (a) The acceptance-corrected \t\ distribution fitted 
with the function f(t) = ae m (solid line), (b), (c), (d) The 
results of the mass-independent PWA (horizontal lines with 
error bars) and a typical mass-dependent fit (solid curve) us- 
ing 0.05 GeV/c 2 mass bins. Only P+ and D+ partial waves 
and their phase difference are shown. The range of the am- 
biguous solutions is plotted with black rectangles, (b) The 
(P+ — D+) phase difference, (c) The intensity distribution 
of the P+ partial wave, (d) The intensity distribution of the 
D+ partial wave. The solid curves in (b), (c), (d) show a 
mass- dependent fit (Fit I) to the P+ and D+ wave intensities 
and the (P+ — D+) phase difference. 



FIG. 3. The results of the mass-independent PWA (hori- 
zontal lines with error bars) and typical mass-dependent fits 
(solid and dashed curves) using 0.1 GeV/c 2 mass bins, (a)-(c) 
show the PWA and Fit 2 (dashed curves) for the P+ - G+ in- 
tensities and their phase difference; (d)-(f) show the PWA 
and Fit 3 (solid curves) for the D+ - G+ intensities and their 
phase difference. The range of the ambiguous solutions is 
plotted with black rectangles. 
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FIG. 4. (a) and (b): Mass-independent PWA for 
0.090 < \t\ < 0.293 GeV 2 /c 2 ; (c) and (d): Mass-independent 
PWA for 0.293 < \t\ < 2.5 GeV 2 /c 2 . The range of the am- 
biguous solutions is plotted with black rectangles. 
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